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ABSTRACT 

Photometric observations made by the NASA Kepler Mission have led to a dramatic increase in 
the number of main-sequence and subgiant stars with detected solar-like oscillations. We present 
an ensemble asteroseismic analysis of 76 solar- type stars. Using frequencies determined from the 
Kepler time-series photometry, we have measured three asteroseismic parameters that characterize 
the oscillations: the large frequency separation (Ai^), the small frequency separation between modes 
of Z = and I = 2 ((5i'o2), and the dimensionless offset (e). These measurements allow us to construct 
asteroseismic diagrams, namely the so-called C-D diagram of 5vq2 versus Az^, and the recently re- 
introduced e diagram. We compare the Kepler results with previously observed solar-type stars and 
with theoretical models. The positions of stars in these diagrams places constraints on their masses 
and ages. Additionally, we confirm the observational relationship between e and T^s that allows for the 
unambiguous determination of radial order and should help resolve the problem of mode identification 
in F stars. 

Subject headings: stars: oscillations 



1. INTRODUCTION 

Solar-like oscillations are global modes that are 
stochastically excited by convection. A revolution in 
their study is now taking place thanks to photometric 
observations bein g made with NASA's Kepler Mission 
(jKoch et al.ll20ldl ). Most of the stars are observed at a 
long cadence of 29.4 minutes, sufficie nt for studying solar - 
like oscillations in red giants (e.g., iHuber et al.ll2Qldl ). 
However, up to 512 stars at a time may be observed 
at a short cadence (SC) of 58.85 s, rapid enough to sam- 
ple oscillations in m ain-sequence and subgiant stars (e.g., 
iChaplin et al.l[20Tot ). This large number of stars makes 
it possible to perform large-scale analyse s in what has 
been t ermed 'ensemble asteroseismology' (jChaplin et al.l 
l2011bD . The aim is to correlate observable global oscil- 
lation parameters with physical properties of the stars, 



such as mass, age and metallicity. This Letter presents a 
first attempt at such an analysis for main-sequence and 
subgiant stars, using a survey of solar-like oscillations 
made during the first ten months of the Kepler Mission. 

2. ASTEROSEISMIC PARAMETERS 

The frequencies of solar-like oscillations for modes of 
high radial order n and low spherical degree I are well ap- 
proxi mated by the so-called asymptotic relation (jTassoull 
[1983): 



5m- 



(1) 



The large frequency separation, Az^, is the spacing be- 
tween modes of the same degree and consecutive order 
and is a probe of the mean stellar density. The small 
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Fig. 1. — Echelle diagram of KIC 1 1395018 showing the frequen- 
cies (black points) as determined bv lMathur et al] 11201 II ). Modes 
with 1 = (circles), 1 = 1 (triangles) and I = 1 (squares) are 
observed. For reference, a gray-scale map showing the power spec- 
trum (smoothed to 1 /.iHz resolution) is included in the background. 
The fits made to the I = and 1 = 2 modes are shown by red lines. 
The values of Ai/ and (ivq2 (at i^niax) and the relationship of e to 
the absolute position of the Z = ridge are indicated by the blue 
arrows, as labeled. 

frequency separations, bv^i^ are sensitive to variations in 
the sound speed gradient near the core in main-sequence 
stars. This is highly dependent upon the mean molecular 
weight of the core and is therefore a probe of the star's 
age. The quantities Ai/ and 8vq2 are labeled in Figure [U 
which shows the ec helle diagram fo r the Kepler subgiant 
KIC 10395018 (.Mathur et al.llMlll . 

The different dependencies of the large and small 
separations have led to the idea of plotting one 
against the other in the so-called C-D diagram 
(jChristensen-Dalsgaardl I1984D . Theoretical calcula- 
tions for the C-D diagram h ave been carried out 
for main-sequence stars (e.e.. iMonteiro et al.l 120021: 
Oti Floranes et all 120051: iMazumdaii 120051: iGai et al.l 
20091 iWhite et al.l 1201 ij) but prior to Kepler, observa- 



tions were only available for a small number of stars. 
An observational C-D diagram has now been produced 
for red giants using Kepl er observations (jBedding et al.l 
l2010at iHuber et al.ll2010( l. In this Letter we present an 
observational C-D diagram for main-sequence and sub- 
giant stars. 

The third asteroseismic parameter in equation ([ij is 
the dimensionless offset e which, as shown in Figure [U 
is related to the absolute position of the I — ridge in 
the echelle diagram. The parameter e has been rather 
neglected un til recently, despite also being investigated 
originally by iChristensen-DalsgaardI ()1984[ ) . An astero- 
seismic diagram based on the relation between e and Av 
was calculated from models by iWhite et al., (,2011. ) from 
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Fig. 2. — A modified H-R Diagram: average large frequency sep- 
aration, Au, against effective temperature for stars in our sample. 
The lines show solar-metallicity ASTEC ( Christcnscn-Dalsgaar^ 
1200811 evolutionary tracks for masses between 0.8 and 1.6 Mq 
in steps of 0.1 Mq. The dashed line indicates the approxi- 
mate location of the c ool edge of the classical instability strip 
llSaio fc GautschvllT998l) 

the zero-age main-sequence to the tip of the red-giant 
branch. They found that evolutionary-model tracks in 
the e diagram are more sensitive to the mass and age of 
subgiants than in the C-D diagram. On the observational 
side, the e dia gram has been publ ished from Kepler red- 
giant data bv IHuber et al.l (|201Cll ). Similar analysis wa s 
performed for CoRoT red giants by lMosser et al.l (|2011| ). 

3. OBSERVATIONS AND DATA ANALYSIS 

The first ten months of Kepler SC asteroseismic obser- 
vations (Q1-Q4) were used for a s urvey, observing abou t 
2000 stars for one month each (jChaplin et al.ll2011bD . 
The time series have bee n prepared from the raw obser- 
vations as described by iJenkins et al.l (l2010f), and fur- 



ther c orrected in the manner described by iGarcia et al.l 
()20lH ). Several methods have been developed to de- 
tect oscillations and extract global oscillation param- 
eters such as Av and z^maxi the frequency of maxi- 



mum power jC ampante et al.l 
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Mosser & Appourcliaux 20 oW- A comparison of these 
methods has been made bv I Verner "eTall (Unni). Not ah 
of the observed solar-type stars showed oscillations in the 
month-long survey data sets, partially due to the impact 
of stellar surface activity on the amplitude of oscillations 
(jChaplin ct al. 2011a). 

After determining which stars showed oscillations and 
their approximate values of Ai^ and fmax, several teams 
estimated the mode frequencies. In general, each star in 
our sample was only fitted by a few teams. The meth- 
ods for extracting frequencies varied betw een teams (see 
ICampante et al.ll2011l : IMathur et aLll201lL and references 
therein for details of the different methods used) . 

Since it is necessary to have clearly identifiable / = 
modes to measure Av and e, and also to have / = 2 
modes to measure 51^02, we are restricted in our sample 
to stars that are bright and/or have the highest ampli- 
tudes. Furthermore, it is necessary to have enough iden- 
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tified modes to make a reliable fit. We required that at 
least three pairs of ^ = and I — 2 modes were identi- 
fied in each star. Stars with ambiguous mode identifi- 
cations were also not considered. These were generally 
hotter stars with large linewidths, that made it difficult 
to distinguish between the ridges of I — 0, 2 and / — 1 
modes (e.g., see I Appourchaux et al.ll20d8l : iBedding et al.l 
l2010bD . 

To select the final frequency list for each star, the 
results from different teams were compared with the 
smoothed power spectrum. Some teams reported modes 
for which there was no peak visible in the power spec- 
trum (usually at low or high frequencies), and those were 
not included. The remaining lists that passed this com- 
parison could provide measurements for Az^, 5vq2 and e 
that we deemed reliable. If more than one frequency list 
remained for a given star, we chose the one that most 
closely fitted the observed peaks in the smoothed power 
spectrum. If no reliable fit could be made, the star was 
excluded from this study. With more Kepler data, which 
will become available for most of these stars, reliable 
measurements of the separations and e should be pos- 
sible in the future. 

Our final sample contains 76 main-sequence and sub- 
giant stars with apparent magnitudes ranging between 
6.9 and 11.7. To illustrate their distribution, we show 
in Figure [2] a modified H-R diagram in which we plot 
(instead of luminosity) against effective temperature. 
Temperatures were determined by An et al. (in prep.) 
from mul ticolor photometry available in the Kepler Input 
Catalog (jBrownet all 1201 ID. Solar-metallic ity ASTEC 
evolutionary tracks (|Christensen-Dalsgaardl I2008D , ne- 
glecting diffusion and core overshoot, are overlaid. Our 
sample spans an approximate mass range of 0.9-1.6 Mq. 
Subgiants dominate the sample because they are intrin- 
sically brighter and their oscillations are higher in am- 
plitude. 

From the frequency list for each star we c alculated At/, 
5vq2 and e using the method described by iWhite et al.l 
([2011). Briefly, a weighted least-squares fit was made to 
the radial (?=0) frequencies to simultaneously determine 
Ai/ and e, as given by equation ([T]). The uncertainties 
in these parameters were determined from the fit. The 
fit was weighted by a Gaussian function centered at the 
frequency of maximum power, fmax- This fitting method 
averages over variation of e with frequency (curvature 
in echelle diagrams) that arise from acoustic glitches 
and t he frequency dependence of the upper turning 
point ([Perez Hernandez fc Christensen-Dalsgaardl 119981 : 
IRoxburghll20'Tol) ! 

In the Sun, 5vq2 is known to decrease linearly with fre- 
quency (Elsworth et al. 1990). In calculating 5vq2 from 
models we performed a linear fit, initially allowing the 
gradient of 5vq2 as a function of frequency to be a free 
parameter. Again, the fit was weighted by a Gaussian 
function and the value of (5:^02 was determined at Vmax- 
However, for the vast majority of stars, the number and 
precision of the frequencies determined from the Kepler 
data do not yet justify the inclusion of the gradient as 
an extra free parameter in the fit. When fitting to the 
data we have kept this gradient constant at the well- 
determined solar value {d5vo2/<^v = —0.0022). In prac- 
tice, keeping this gradient constant did not significantly 
affect the measured value of (5 1^02, or its uncertainty. 



4. RESULTS 
4.1. C-D diagram 

Figure [3| shows the C-D diagram for our sample 
(filled red circles). We also include: (i) five Kepler 
stars with previously published frequency lists (open 
red circles) - KIC 11026764 (known within the Ke- 
pler Asterosei s mic S cience Consortium as 'Gemma'; 
IMetcalfe et all IMol) . KIC 10273246 ('Mulder') and 
KIC 10920273 ('Scully'; both lCampante et aLllMl . and 
KIC 11234888 ('Tigger') and KIC 11395018 ('Boogie'; 
both lMathur et al.l l2011l: (u ) three CoRoT s tars (green 
triangles) - HD 49385 (Dche uvels eFalllmol). HEi 4993 3 
([Benomar et al.|[200l ) and HD 52265 ([Ballot et al.|[20Tll) : 
(iii) six targets obser ved from the ground (black di- 
amonds) - a Cen A ([Bedding et al.ri2004D, a Cen B 
([Kjeldsen et al.' '2005'), /3 Hyi ([Bedding et al.' '20071), 
77 Boo (Kieldscn et al. 2003), u ArafBouchv ct al. 2005|) 
and r Get (Tcixcira ct al. 2009); (iv) the Sun (usual sym- 
bol; iBroomhjiJljjtj^^ and (v) 470 Kepler red giants 
for which Ai/, 5yQ2 and e were measured by Huber et al. 
(2010; black circles). 

Stars evolve in the C-D diagram from the top-right to 
the bottom- left. In Figure [3] we al so show t he model evo- 
lutionary tracks as determined bv lWhite~e t al. (2011) for 
near-solar metallicity {Zq = 0.017), which are based on 
the ASTEC stellar models shown in Figure [H During 
the main-sequence, the evolutionary tracks are well sep- 
arated (at least for masses below ~ I.5M0) but converge 
during the subgiant and red giant phases. Decreasing the 
metallicity has the effect of shifting the tracks up and to 
the left, and in the opposite direction for higher metallic- 
ity. We show the tracks for metal-poor {Zq = 0.014) and 
metal-ric h (Zo = 0.022) so lar-mass models for compar- 
ison (see IWhite et~all[20nl for model C-D diagrams of 
different metallicities) . To ensure that the observed val- 
ues are directly comparable to those derived from models, 
the same recipe for measuring Av, e and 5i^o2 from data 
as described in Section [3] was also applied to the models. 

The positions of the Kepler stars in the diagram are 
well covered by the models. In particular, the data show 
the convergence of different masses during the subgiant 
phase, as predicted by the modeling. As a result of 
this convergence, the value of Svo2 must be measured 
with high precision for subgiant and red giant stars if 
it is to provide significantly more information about the 
structure of the star than Av alone. This was noted in 
the modeling of the Kepler subgiant KIC 11026764 by 
IMetcalfe etall (|2010( ). To better constrain the masses 
and ages of subgiants, the e diagram may prove more 
useful (see Section 1421 below) . 

The presence of Z = 2 mixed modes in the oscilla- 
tion spectrum can complicate the measurement of 5vq2 
by shifting modes from their asymptotically expected fre- 
quencies. This is caused by coupling betw een jp modes 
in the envelope and q modes in the core ([(^sakil[T975l : 
lAizenman et al.lll977( ). The coupling for Z = 2 is con- 
siderably weaker than for I = 1, but can still result in 
a significant shift in the measured value of 5vq2- This is 
one cause of the scatter for subgiants of the observed val- 
ues and also of the model tracks, which further increases 
the difficulty of determining stellar properties from the 
values of Av and 5vq2 alone. Metallicity also contributes 
to the observed scatter. 
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Fig. 3. — C-D diagram, showing 5uq2 versus A.u. The filled red circles are the 76 Kepler stars from this work, while open red circl es are 
Kepler stars with previously published frequency lists (see text). Also shown are Kepler red giants (black circles: [Huber et al.ll20l"0l '). and 
main-sequence and subgiant stars from CoRoT (green triangles; see text) and ground- based observations (black diamonds; see text). The 
Sun is marked by its usual symbol. Error bars show the uncertainties derived from the linear fits for both h.v and 51^02, although those in 
At' are generally too small to be visible. Model tracks for near-solar mctallicity (Zq = 0.017) increase in mass by 0.1 Mq from 0.8 Mq to 
2.0 Mq (light blue to red lines). Also shown are tracks for metal-poor {Zq = 0.014; dotted) and metal-rich {Zq = 0.022; dashed) solar-mass 
models. The section of the evolutionary tracks in which the models have a higher Teg than the approximate cool edge of the classical 
instability strip ( Saio & Gautschy 1998;) are gray: they are not expected to show solar-like oscillations. Dotted black lines are isochrones, 
increasing from Gyr (ZAMS) at the top to 13 Gyr at the bottom. 



The C-D diagram is primarily useful for main-sequence 
stars, where the evolutionary tracks are well-separated. 
If the position of a star in the C-D diagram can be well 
determined and the metallicity is accurately estimated, 
then its mass and age can be constrained. The physical 
process es incorporated in to the models can also affect the 
tracks (|Mazumdaii I2005D and so masses and ages deter- 
mined from the C-D diagram will be model dependent. 
Measurements of temperature and lumin osity may aid in 
distin guishing between different models (jMonteiro et al.l 
[2001 . 

For the stars in our sample, we do not yet have reli- 
able metallicities, so we are unable to provide accurate 
determinations of their masses and ages based upon their 
position in the C-D diagram. However, we are able to 
consider the precision with which these quantities could 
be determined if we assume a particular uncertainty in 
metallicity. We have considered several main-sequence 
stars of varying mass and age. With an uncertainty in 
[Fe/H] of ±0.1, the masses of these stars could be de- 
termined to within 4 to 7%. Age can generally be de- 
termined to within 1 Gyr. The precision of the age is 
heavily influenced by the uncertainty in 5vq2 and the po- 
sition of the star in the C-D diagram. Many of these stars 



are being observed for an extended period of time with 
Kepler and the precision of these measurements will un- 
doubtedly improve. With supporting ground-based spec- 
troscopic observations and detailed modeling involving 
the individual frequencies, the fundamental properties 
of these stars will become well-determined and the data 
should provide significant tests of stellar models. 

4.2. The e Diagram 

Figure m shows the observational e diagram. Observa- 
tions cannot determine the radial order n directly, and 
so it is possible for e to be uncertain by ±1, particu- 
larly in the subgiants {/S.v in the range 20-80 /.tHz). For 
these, we have taken e to be in the range 0.7-1.7, but 
note there is some ambiguity for stars near the extremes 
of this range. The measurement of e is complicated by 
its close relationship to ls.v: a small change in Aj/ can 
induce a large change in e. We have also measured e 
using an alternative method: the variation of the large 
separation wi t h freq uency was measured as described by 
iMosser et al.l (|2010[ ). before the radial modes were glob- 
ally fit and e derived, taking into account the mean cur- 
vature. A comparison of the values obtained by the two 
methods showed good agreement, although small system- 
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Fig. 4. — The e diagram, with near-solar metallicity {Zq = 0.017) model tracks. Tracks increase in mass by 0.1 Mq from 0.8 to 
2.0 Mq. Also shown are tracks for metal-poor (Zq = 0.014; dotted) and metal-rich (Zq = 0.022; dashed) solar-mass models. All colors, 
lines and symbols are the same as for Figure O although, for clarity, models in the instability strip are not shown. Error bars show the 
uncertainties from the simultaneous linear fitting of /\v and e (the uncertainties in l\v are mostly too small to see). 



atic offsets exist, typically about 0.1. This offset is prob- 
ably due to the combined effects of curvature (departure 
from equation ([1])) and the slightly different range of fre- 
quencies over which /S.v and e were measured. A single 
method must be applied to both models and data used 
to ensure consistency. In this Letter we have used the 
method outlined in Section [3] 

The observed stars in Figure |4] are offset to the right 
of the models. This offset is well-known from helioseis- 
mology, in which there is a discrepancy between the ob- 
served and computed oscillation frequencies of the Sun 
arising from improper modeling of the near-surface layers 



(iDziembowski et all 119881 : iChristensen-Dalsgaard et all 
I1996D . A rigorous comparison of the observations with 
models requires that the offset be taken into account, ei- 
ther by a proper modeling of near-surface layers, or by 
applying an empirical correctio n to the models as sug- 
gested bv lKieldsen et al.l ()2008[ ). It remains to be shown 
whether this empirical correction is valid for stars over a 
wide range of evo lutionary sta. t es. 

As discussed bv I White et al.l (|2011[ ). the value of e may 
be a significant asteroseismic constraint on the masses 
and ages of stars. This is particularly useful for sub- 
giants, for which the convergence of the evolutionary 
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Fig. 5. — e as a function of effective temperature, witfi near-solar 
metallicity (Zq = 0.017) model tracks. Tracks increase in mass 
by 0.1 Mq from 0.8 to 2.0 Mq. All colors, lines and symbols 
are the same as for Figure [S] although, for clarity, models in the 
instability strip are not shown. 

tracks in the C-D diagram limits the diagnostic potential 
of the IS.v-5vq2 relation. However, it can be difficult in 
practice to constrain e b ecause it var i es wit h frequency, 
as discussed in detail by I White et al.l (|2011f) . As appar- 
ent in Figure |31 the value of e is sometimes poorly deter- 
mined. It should also be s tressed that the m odel tracks 
will vary with metallicity (|White et al.ll20H"l ). as can be 
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seen from the solar-mass tracks of different metallicities 
in Figure |4l Nevertheless, for many stars e is well con- 
strained and with supporting spectroscopy to constrain 
metallicity, the e diagram shows potential as a useful as- 

teroseismic diagnostic tool. 

Finally we note that IWhite et alj (|2011[ ) demonstrated 
a relationship between e and effective temperature. We 
confirm this relationship observationally with our larger 
sample (Figure [SJ. With this diagram, the ±1 ambiguity 
over the value of e, and therefore radial order, can be 
resolved. The extension of this relationship to higher 
Tcff shows promise for resolving the mode identification 
problem in F stars referred to in Section [3] (White et al. 
in prep). 
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